Phonon Softening and Elastic Instabilities 
in the Cubic-to-Orthorhombic Structural Transition of CsH 
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The cubic-to-orthorhombic structural transition occurring in CsH at a pressure of about 17 GPa 
is studied by ab initio calculations. The relative stability of the competing structures and the 
transition pressure are correctly predicted. We show that this pressure-induced first-order transition 
is intimately related to a displacive second-order transition which would occur upon application of 
a shear strain to the (110) planes. The resulting instability is rationalized in terms of the pressure- 
induced modifications of the electronic structure. 
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In recent years the availability of the diamond anvil- 
cell technology has renewed the interest in the high- 
pressure properties of alkali halides which are the sim- 
plest and prototypical among ionic solids. In the case of 
cesium halides, these studies have been carried out espe- 
cially in the search of band-overlap metallization which 
occurs at pressures in the Mbar range. In the quest 
of metallization, unexpected phase transformations to 
tetragonal and orthorhombic structures have been ob- 
served experimentally |^;^,^,|[|^ and studied theoreti- 
cally 

Metal hydrides are considerably more covalent than 
the corresponding halides. In fact, transition-metal hy- 
drides are considered to be covalent materials, whereas 
alkali hydrides are structurally rather similar to the cor- 
responding halides. Understanding the behavior of these 
materials at high pressure would presumably allow to 
better characterize the bonding properties of hydrogen 
in extreme pressure conditions and thus provide valuable 
hints towards the long searched goal of hydrogen metal- 
lization. 

Alkali hydrides at zero pressure crystallize in the rock- 
salt (cubic Bl) structure, while most of them undergo 
a transition to the cubic B2 (CsCl-likc) structure at an 
applied pressure of a few GPa [ p^pl|jl^ . Recently, a 
second transition from the B2 structure to a new or- 
thorhombic phase has been observed to occur in CsH at 
an applied pressure of about 17 GPa jl^. The new phase 
has been assigned the CrB structure which belongs to the 
space group. 

In order to assess the driving mechanisms of the tran- 
sition we have undertaken a series of first-principles cal- 
culations of the electronic, structural, elastic, and vibra- 
tional properties of CsH as well as of their dependence 
upon pressure. To this end, we have employed density- 
functional-theory (DFT) within the local-density approx- 
imation (LDA) pj], and its gradient-corrected (GC) gen- 
eralizations Jist . Our calculations have been performed 
using norm-conserving pseudopotentials p6[ . In the case 



TABLE I. Comparison between calculated and experimen- 
tally observed structural properties of CsH. ao (a.u.) is the 
equilibrium lattice parameter in the Bl phase, Bo (GPa) the 
bulk modulus. Pi* and P2 the transition pressures (GPa) to 
the Bl and CrB structures respectively, while AVi and AV2 
are the corresponding volume changes. Vy/Vo is the ratio of 
the volume at the B2 —* CrB transition to the equilibrium 
volume. 
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of Cs, 5s and bp semi-core states have been treated as va- 
lence states. For LDA calculations we have used the same 
Cs potential as in Ref . , while Hydrogen is described by 
a norm-conserving pseudopotential constructed for the Is 
wave- function, so as to smoothen it somewhat and make 
our calculations well converged with the 25 Ry kinetic- 
energy cutoff which we use. New potentials have been 
generated for GC calculations. The vibrational proper- 
ties are determined using the density-functional pertur- 
bation theory (DFPT) described in Ref. 0. 

In Table | we report our results for some structural 
properties of CsH in the three phases studied here: Bl, 
B2, and CrB. The LDA fails to account for the stabil- 
ity of the rocksalt structure at zero pressure (negative 
Bl B2 transition pressure), while GC calculations 
predict the correct sequence of transitions. GC-DFT 
also considerably improves the agreement between cal- 
culated and observed equilibrium lattice parameters and 
bulk modidi. The inclusion of Cs semi-core states in the 
valence manifold implies that the pseudopotential trans- 
ferability is optimal around the corresponding orbital en- 
ergies, rather than at the energies of the valence states. 
As a consequence, this treatment of semi-core states — 
although essential to obtain sensible results ^ — might 
require the use of more than one reference state in order 



1 




0.0 



4.0 



8.0 



a 

12.0 



15.8 



19.5 



FIG. 1. Decomposition of the distortion leading to the CrB 
structure (lower panel) into an elastic (upper left) and a vi- 
brational (upper right) contributions. Large and small circles 
represent Cs and H atoms respectively. Dark and light shad- 
ing refer to the topmost and second (110) layers respectively 



to ensure optimum transferability. Whether or not the 
relatively poor quality of the LDA predictions is due to 
deficiencies in the Cs pseudopotential and the improve- 
ments achieved by using GC-DFT due to a fortuitous 
cancellation of errors is a matter which deserves further 
investigations. In any events, our results seem to confirm 
the importance of the gradient corrections to the LDA in 
the description of structural phase transitions, recently 
claimed in the case of the diamond to /3-tin transition in 
silicon Both LDA and GC-DFT calculations cor- 

rectly predict that a transition from the B2 to the CrB 
phases would occur at pressures somewhat larger than 
10 GPa. In this case too, GC-DFT gives a value of the 
transition pressure in closer agreement with experiments 
(15 vs. 17 GPa). Note however that — due to hystere- 
sis effects — the experimental estimate done when loading 
the sample only provides an upper limit to the tran- 
sition pressure. Also the volume discontinuity and the 
ratio between the transition volume and the equilibrium 
one are well reproduced by our calculations. 

The CrB structure can be viewed as resulting from 
a continuous deformation of the cubic B2 structure in 
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FIG. 2. Enthalpy vs. angle of deformation for different val- 
ues of the applied external pressure. In the inset we show the 
correponding values of the amplitude of the lattice distortion, 
u. 



which a shear of the cubic edges in the (110) planes is 
coupled to a vibrational distortion of the atomic lattice, 
having the periodicity of a zone-border transverse phonon 
at the M point of the Brillouin zone (BZ), q = 5, 0), 
with symmetry M^, and polarized along (001). In Fig. 
1 we display how the distorted (110) plane (lower panel) 
results from the combination of the elastic (upper left) 
and the vibrational (upper right) deformations. In the 
ideal CrB structure the a and h crystallographic parame- 
ters indicated in Fig. 1 are in the ratio b/a — 2\f2 which 
would be appropriate to a cubic structure. Due to the 
lower symmetry of the CrB structure with respect to the 
cubic one, the actual value of hja slightly differs from 
the ideal one and depends somewhat on pressure. This 
dependence is however very weak: at the transition pres- 
sure, for instance, one has hja k, 2.9. The third crystal- 
lographic axis is orthogonal to the ah plane, and the cor- 
responding crystallographic parameter has an ideal value 
c = ^/2a. The pressure dependence of c/a is somewhat 
stronger than that of hja: at the transition, for instance, 
one has c/a « 1.3. 

In Fig. 2 we show the T = enthalpy, H = E + PV, 
of the system constrained to a given elastic deformation, 
as a function of its amplitude which we measure through 
the parameter t — tan(a) x b/a [|l^. We see that at any 
positive pressure a second minimum exists for t — 1, cor- 
responding to the CrB orthorhombic structure. Actually, 
the exact value of a slightly depends on pressure because 
so does b/a. For pressures higher than P2 ~ 11 GPa 
the second minimum becomes more stable giving rise to 
a first-order phase transition. The amplitude of the vi- 
brational distortion, u = x/a, which minimizes the en- 
thalpy, Umin, depends on t, but for f = 1 it is constrained 
by symmetry to be Umin = 1/4, independent of pressure. 
The dependence of Umin upon the amplitude of the elas- 
tic distortion, t, is shown in the inset of Fig. 2. We see 
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FIG. 3. Frequency of the phonon mode (cm ^) as a function of the apphed hydrostatic pressure for different angles of 
deformation of the cubic structure (left), and as a function of the angle of deformation for different pressures (right). 

u direction, while nothing can be concluded in principle 



that Umin remains equal to zero up to a critical value of 
t ~ t* which depends on pressure, and that it saturates 
to 1/4 for t = 1. This behavior is typical of the order 
parameter at a second-order phase transition, and it in- 
dicates therefore that a phonon frequency softens when 
the elastic distortion becomes larger than t* . In Fig. 3 
we display the M2 phonon frequency as a function of t at 
different pressures, P, and as a function of P for different 
elastic distortions. 

The softening of the M2 frequency in correspondence 
to some value of the elastic distortion could have been ex- 
pected on the basis of symmetry considerations. To see 
this, let us consider the crystal energy as a function of t 
and u, E{t, u). Crystal symmetry requires that E(t, u) = 
E{±t, u) = E{t, ±u) = E{t + 2,u+\) = E{t + 2,u-\). 
These relations imply that i?(t, u) is stationary at the 
point (i, m) = (2,0). This point lies in between the four 
equivalent minima at (0,0), (2,±i), and (4,0). Assum- 
ing that the energy landscape is simple, we arrive at the 
conclusion that (2, 0) is a maximum. In particular, one 
has that d'^Ejd'u} < 0, and hence there exists a value 
of < t < 2 for which ul - d'^E/du^ = 0. The or- 

thorhombic phase corresponds to the (1, j) point in the 
(i, u) plane which is also stationary because of symme- 
try. This point is located halfway between the two points 
(1,0) and (1, ^) which are equivalent by symmetry and 
stationary with respect to variations of u {dE/du = 0). 

If t* were larger than 1, d^E/du^ would be positive at 
(1, 0) and (1, i), and therefore the assumption of a simple 
energy landscape would imply that (1, ^) is a maximum 
with respect to u, and the orthorhombic structure unsta- 
ble. It is easy to see that a first-order transition with a 
discontinuous variation of u would occur in this case as 
the elastic distortion crosses t ~ 1. In fact, in this case 
the line {t, 0) is stable with respect to variations of u, 
up to t = 1. For t > 1 one has that E(t, i) < E{t,0). 
Hence, a discontinuous jump in u would occur a,t t — 1 
where Eit, i) = E{t,0). 

Arguments similar to those expounded above show 
that for t* < 1 the (1, j) point is a minimum in the 



for other directions, but the fact that local stability is 
allowed. Our calculations indicate that in the present 
case t* < 1 and that the CrB structure is indeed locally 
stable. For pressures higher than 11 GPa (15 GPa within 
GC), this structure is favored with respect to the cubic 
one, giving thus rise to a first-order transition. 

The mechanisms driving the softening of the M2 
phonon are related to the metallization of CsH under an 
applied pressure and/or shear deformation. At equilib- 
rium {P = 0, t = 0), the maximum of the valence band 
is practically degenerate between the R = (555) and 
X = (OOi) points of the BZ, and the minimum of the 
conduction bands lies at the R point. A small applied 
pressure lifts the degeneracy between the valence-band 
maxima, rising the X point with respect to R. The re- 
sulting fundamental band gap, Eg, corresponds therefore 
to a X R transition, whose transferred momentum is 
M = (^^0) , i.e. the wave- vector of the phonon which 
goes soft. 

In Fig. 4 we display the electronic energy bands of CsH 
at P = 17.5 GPa, for t = and for t = 0.3. In the dis- 
torted structure — corresponding to a non- vanishing value 
of t — the R and X points are folded to a same point of 
the BZ, and so are the F and M points. The gap, which 
is indirect in cubic symmetry, becomes thus direct in the 
distorted structure. We see that the elastic distortion 
considerably favors metallization, and thus enhances the 
screening of the phonon frequencies at the M point of 
the cubic BZ. 

In Fig. 5 we display the dependence of Eg upon pres- 
sure and deformation angle. At small pressures, the 
phonon frequencies softens almost for the same values of 
t for which the gap closes, whereas for pressure larger 
than « 18 GPa, it takes a larger value of t to soften the 
phonon frequency than to close the gap. If the bands 
were parabolic around the valence- and conduction-band 
edges, the phonon would go soft at the metallization. 
In fact, the independent-electron polarizability, Xo(QAf) 
would diverge logarithmically in this case ||2^], and the 
restoring force for a lattice distortion of wave- vector 
would correspondingly vanish. If the valence and con- 
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FIG. 4. Electron energy bands of CsH close the fundamental energy gap (which is evidentiated by the shadowed area) at an 
applied pressure of 17.5 GPa in the cubic structure (left) and under a finite shear deformation (right) 
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FIG. 5. Fundamental energy gap (eV) as a function of the applied hydrostatic pressure for different angles of deformation of 
the cubic structure (left), and as a function of the angle of deformation for different pressures (right). 



duction bands have different shapes, the logarithmic di- 
vergence would no longer hold, and all that can be said 
in this case is that the more these shapes are similar, the 
more metallization contributes to the phonon softening. 
In any events, when the amplitude of the elastic distor- 
tion is larger than t*, the spontaneous vibrational distor- 
tion due to phonon softening re-opens the band gap and 
makes the orthorhombic phase insulating, in agreement 
with the behavior of Eg found in 
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